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ABSTRACT 

I study the behaviour of the maximum rms fractional amplitude, r max and the max- 
imum coherence, (2 max , of the kilohertz quasi-periodic oscillations (kHz QPOs) in a 
dozen low-mass X-ray binaries. I find that: (i) The maximum rms amplitudes of the 
lower and the upper kHz QPO, r^ ax and r^ &x , respectively, decrease more or less 
exponentially with increasing luminosity of the source; (ii) the maximum coherence 
of the lower kHz QPO, Q max , first increases and then decreases exponentially with 
luminosity, at a faster rate than both r^ ax and r^^; (iii) the maximum coherence 
of the upper kHz QPO, Q max , is more or less independent of luminosity; and (iv) 
'"max and Q max show the opposite behaviour with hardness of the source, consistent 
with the fact that there is a general anticorrelation between luminosity and spectral 
hardness in these sources. Both r max and Q m ax in the sample of sources, and the rms 
amplitude and coherence of the kHz QPOs in individual sources show a similar be- 
haviour with hardness. This similarity argues against the interpretation that the drop 
of coherence and rms amplitude of the lower kHz QPO at high QPO frequencies in 
individual sources is a signature of the innermost stable circular orbit around a neu- 
tron star. I discuss possible interpretations of these results in terms of the modulation 
mechanisms that may be responsible for the observed variability. 
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1 INTRODUCTION 

Kilohertz quasi-periodic oscillations (kHz QPOs) in the X- 
ray flux of low-mass X-ray binaries have drawn much atten- 
tion since their discovery, about ten years ago. The reason 
for this continued interest is that since they most likely re- 
flect the motion of matter very close to the neutron star 
(or black hole) primary in these systems, these QPOs may 
provide one of the few direct ways of measuring effects that 
are unique to the strong gravitational-field regime in these 
systems. In general two simultaneous kHz QPOs are seen 
in the power density spectra of low-mass X-ray binaries, the 
lower and the upper kHz QPO according to how they appear 
sorted in frequency. 

Most ( but not all) of the work on QPOs in these 
years (see Ivan der Klisl l2006t for a review) has fo- 
cused on the frequencies of these QPOs, because those 
frequencies provide insights into the dynamics of the 
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system. Several models have been proposed to ex- 
plain the observed relation between both kHz QPO 
frequencies jMiller. Lamb, fc Psaltid Il998t lLamb fc Millei 



20011: IStella fc Vietril ll99SL Il999l: lOsherovich fc Titarchuk 



119991: iTitarchuk fc Osherovicbf 1999; lAbramowicz et al 
booa iRebuscol l2005t , their relation wi t h the spi n 



of 



the n eutron star (e.g., iMiller et alJ Il998t iKluzniak et alJ 
2004) as well as with other low- frequency variabil- 



hy (e.g., 


IStella fc VietrJ 


1998; 


Osherovich fc Titarchuk 


1999: [Titarchuk & Osherovichlll999l: ITitarchuk et al.ll200ll: 



broad-band variability. 

From the beginning there has been interest on 
the other two properties of th e kHz QPOs, the i r am- 
plitude and coherenc e (e.g., Ivan der Klis et all Il997l : 
IWiinands et alJ Il998afl : but systematic studies of those 
other QPO prope rt ies st ar ted t o take off slightl y later 
on Jjonker et al.l EM 120021: iDi Salvo et all l200l|: 
iDi Salvo. Mendez. fc van der Klisl 120031 1 van Straaten et all 
hood, l2002t iMendez. van der Klis. fc Fordl 1200 ll : 
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Figure 1. The coherence of the lower (filled black symbols) and 
upper (open gray symbols) kHz QPO as a function of the QPO 
frequency for the sources studied in this paper. Notice that the 
range spanned by the y-axis decreases from the panels at the top 
to those at the bottom. 
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Figure 2. The rms amplitude of the lower (filled black symbols) 
and upper (open gray symbols) kHz QPO as a function of the 
QPO frequency for the sources studied in this paper. To avoid 
having to compress the scale of the y-axis for some sources, the 
ordering of sources in this Figure is different from that in FigurelTl 
There are no measurements of the true rms amplitude of the 
QPOs in Sco X-l (see § 12.41 . and therefore I do not plot this 
source in this Figure. 



van Straaten. van der Klis. fc Mended 12 0031: iH oman et alJ 
20021: iKuznetsovl 1200* IfWt et al.l l2005al: 

Barret. Olive, fc Miller 2005b d lAltamirano et alJl200f3) . 

Recently. iBarret et alJ i2005allbl3) carried out a system- 
atic study of the kHz QPO coherence and rms amplitude in 
three X-ray binaries, 4U 1636-53, 4U 1608-52, and 4U 1735- 
44. They find that in all three sources the coherence and rms 
amplitude of the lower kHz QPO increase slowly with fre- 
quency, and after they reach their maximum values they 
decrease abruptly as the QPO frequency keeps on increas- 
ing. (The sudden drop is most noticeable in the coherence 
of the lower kHz QPO; in the case of the rms amplitude the 
drop is less abrupt.) They propose that this behaviour is 
due to effects related to the innermost stable circular orbit, 
ISCO, around the neutron star in these systems. 

Triggered by these results, in this paper I investigate 
the dependence of the maximum coherence and rms ampli- 
tude of both kHz QPOs in a large sample of low-mass X-ray 
binaries. In J5| I describe the data I use in the rest of the 
paper. Since I collected most of the data from the literature, 
I spend some time describing how those data were obtained 
and the selection criteria. In Sj3]l present the results, and I 
discuss them in <0] 



2 DATA SELECTION AND ANALYSIS 
2.1 Data selection 

All the data that I use in this paper were obtained over the 
last 10 years with t he Proportional Counting Array, (PCA; 
Ijahoda et alJll996l). on board the Rossi X-ray T iming Ex- 
plorer (flXrg; lBradt. Rothschild, fc Swanklll993l) . 

I collected most of these data from the literature. For 
this I searched, for as many sources as possible, all published 
values of the rms fractional amplitude and the coherence 
of the kHz QPOs. The coherence Q of a QPO is defined 
as Q — fQpo/X, where ^qpo an d ^ are the frequency and 
the full- width at half-maximum of the QPO. Some authors 
report A and v instead of Q; in those cases I calculate Q using 
the previous formula. The rms amplitude, r, is calculated 
from P, the integral from to 00 of the Fourier power u nder 
the QPO, and the source intensity, S, as r = 100 x -J P/S 
{van der Klisl H9951: from this definition, r is expressed as a 
percent of the total intensity of the source. 

From now on I will us e the naming convention of the 
kHz Q POs introduced by iBelloni. Psaltis. fc van der Klisl 
(2002), in which the lower kHz QPO is called Li, and the 
upper kHz QPO is called L u . The frequency, coherence and 
amplitude of each QPO carry a subscript or superscript I or 
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u, respectively. A clarification is needed abou t what I call 
kHz QPO in this paper. Some authors (e.g., IPsaltis et alJ 
1999) identify low-frequency features as possibly being Le 
on the basis of an extension of the QPO frequency-frequency 
correlation, v n — vt, to low frequencies (below vt ~ 50 Hz). It 
remains still uncertain whether these id entifications are cor- 
rect fcf lvan Straaten et al]l2003t 12005). In this paper I use 
the expression "kilohertz QPO" to refer to features in the 
power spectra of neutron star systems that have frequencies 
> 150 Hz, and that have not been identified as hectohertz 
QPO fe.g- lvan Straaten et al]l2000l) . 

Since the goal is to compare the maximum values of Q 
and r of each kHz QPO in different sources, one of the re- 
quirements in collecting the data is that there are sufficient 
measurements of both parameters to have confidence that 
the maximum values have been measured. This is in part 
facilitated by the fact that we know the general dependence 
of Q and r on QPO frequency for both QPOs. For instance, 
van Straaten et al Jll2000l.l2q02 j .l2003fl.lDi Salvo et alJfcOOll 



and 



20031). IMendez et alJ dpOOlf ). Barret et alJ ( 12005 
Alt^mfiran^^t^dTl200a) show that both n and Qe increase 
with Vt, they peak at intermediate to high values of ve, and 
then they decrease as ve increases ( see e.g., Figure 3 i n 
iDi Salvo et all l200ll and Figure 3 in iBarret et"afl l2005bl) . 
On the other hand, r u is more or less flat at low values of 
v u where it peaks, and then decreases as v u increases, while 
Q u is more or less constant with v n , although it increases to- 
ward s higher u u value s (see e.g., Figure 3 in [ Pi Salvo et ail 
l200ll. and Figure 10 in lAltamirano et al.l2006l) . For this rea- 
son, here I compile data from papers that present systematic 
studies of the QPO properties, or that have enough measure- 
ments to deduce the systematic trends of those properties. I 
do not use result from papers reporting the discovery of kHz 
QPOs in a source, in which only one or two measurements 
of the QPO properties are available. 

In Figure Q and Figure |2] I present the plots of Q 
and r vs. QPO frequency for all the sources described in 
this paper: 4U 1608-5 2 jBarret et al. l2005al). 4U 1636 - 



52 JBarret et al.ll2005bD 4U 1 728-34 iBarret et ai]l2006l) 
Aql X-l Wendez et ahl I2001T) 4U 1820-3 barret et alJ 
20061 ). 4U 1735-3 4 iBarret etHI l2006h. 4U 0614+09 




and GX 340+0 iJonker et aDExwh . As it is apparent from 
this figure, in most cases Qe, re, and r u have well-defined 
maxima that occur at QPO frequencies in between the min- 
imum and maximum frequency observed for each QPO in 
each source. For three sources (GX 17+2, GX 5-1, and GX 
340+0), however, the maximum value of Qe occurs at the 
edge of the frequency range. Similarly, for two sources (GX 
5-1 and GX 340+0), the maximum value of re, and for three 
sources (4U 1820-30, Cyg X-2, and GX 340+0) the maxi- 
mum value of r u occurs at the edge of the frequency range. 
For Qu, it is more often the case that the maximum ob- 
served value occurs at the edge of the frequency range (five 
out of eleven sources). While it may still be possible that in 
these cases the maximum value of Q or r has not yet been 
observed, on the basis of the slopes of the relations of Q and 
r with v in Figure Q and it seems very unlikely that the 
(unseen) maxima are too different from the largest values 
so far measured. Except the case of re in GX 5-1, for which 



the maximum value might still be somewhat larger than the 
largest value so far observed, for the other sources Q and r 
appear to level off at the edge of the frequency range. Based 
on this, I estimate that in those cases the maximum value 
can be at most ~ 20% higher than the value I use in this pa- 
per; this potential difference has no effect on the conclusions 
of the paper. Notice also that since the signific ance at which 
a QP O is detected is proportional to r 2 x Q 1 ^ 2 ivan der Klisl 
1989), if Q or r were significantly larger than so far observed, 
QPOs would have most likely been detected above (below) 
the current maximum (minimum) frequencies. 

A word must be said about the different ways in which 
Q and r are measured in these sources, and how these dif- 
ferent ways of measuring them can affect their values: 

(i) The first QPO measurements were obtained over con- 
tinuous RXTE observations, usually stretches of up to ~ 
3000-s exposure time interrupted by the occultation of the 
source below the Earth's horizon, or by the passage of the 
satellite by the South Atlantic Anomaly, when many in- 
struments were switched off for safety reasons. Because the 
QPO frequency typic ally changes by a fe w tens of Hz within 
a thousand seconds iBerger et al.lll99^) . the QPOs appear 
broader than they actually are in the average power spectra 
of those observations. Since the width measured in the av- 
erage power spectrum, A b s , when a QPO of intrinsic width 
A drifts by an amount Sv is roughly A obs ~ ^\ 2 + {Sis) 2 , 
clearly intrinsically narrower peaks are more affected by the 
QPO frequency drift than intrinsically broader peaks. Be- 
cause in general Le is narrower than L u , the width of Le is 
the most affected. On the other hand, since the total power 
in the QPO is conserved in the averaging process, the rms 
amplitudes of the QPOs are less affected by the QPO fre- 
quency drift. (Actually, the drift changes the shape of the 
QPO in the average power spectrum, and because the func- 
tion used to fit them, usually a Lorentzian, no longer rep- 
resents them properly, also the rms amplitudes are affected, 
although to a lesser degree) . 

(ii) When more data per source were collected, and it 
was realised that the QPO frequencies depended on other 
source parameters, it became clear that a much better way 
of studying the QPO properties was to average data taken 
at different epochs but having similar range of values of 
those other source parameters. The typical examples are 
data selection based on the colours of the source, or more 
specifically the position along the track traced out by the 
source in a colour-colour or colour-intensity diagram (e.g. 
Jonkere^il.l [2QQ^|: [Horn^n e^jj 1200^ Ivan^S^raa^eji^^d 
| 200(1 l200alDi Salvo et alJ l200ll 120031: lAltamirano etaT 
l2005l 120061) . or the frequency of one of the kHz QPO usu- 
ally Le fe.g.. IMendez et al . I I2OO1I IBarret et alJl2005bl ). but 



endez 



sometimes also L u (Me ndez fc van der Klis 2000) , or the fre - 
quency of a low-freauencv QPO (e.g.. Ijonker et alj [2000). 
While these selections tend to reduce the spurious broad- 
ening of the QPOs due to the QPO frequency drift, this 
effect may still be present in the results, especially when the 
parameter used in selection spans a large interval. 

In this paper I try to use as much as possible data col- 
lected according to the second method just described, es- 
pecially for the values of Q, and in particular for Qe, which 
would be the quantity that is most affected by the averaging. 

In a few cases one or more of the QPO parameters 
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Table 1. Maximum coherence and maximum rms amplitude of the lower kHz QPO for the sources studied in 
this paper. 



Source 



(Hz) 



(L/L Edd )q^ 



(%) 



(Hz) 



(L/ X Edd )_{ 



Ref 



4U 0614+09 


28.5+ 4.0 


730 


0.065 + 0.012 


9.3 ± 0.3 


610 


0.065 + 0.016 


1,2,3,4 


4U 1608-52 


247.0 + 16.0 


830 


0.040 + 0.010 


10.2 ± 0.4 


670 


0.030 ± 0.075 


4,5,6,7,8 


4U 1636-53 


248.0 + 17.0 


840 


0.085 + 0.021 


8.9 ± 0.2 


780 


0.085 + 0.021 


4,9,10,11 


4U 1728-34 


188.0 + 18.0 


900 


0.080 ± 0.020 


9.2 ± 0.6 


760 


0.070 + 0.017 


4,5,7,12 


4U 1735-44 


130.0 ± 12.0 


770 


0.120 + 0.030 


7.3 + 0.2 


710 


0.120 + 0.030 


4,13,14 


4U 1820-30 


110.0 + 17.0 


700 


0.200 + 0.050 


5.9 + 0.3 


600 


0.200 ± 0.050 


4,5 


Aql X-l 


190.0 + 19.0 


850 


0.019 + 0.047 


8.7 + 0.6 


760 


0.015 + 0.041 


4,5,7 


GX 5-1 


13.7+ 4.9 


620 


1.800 ± 0.450 


3.6 + 0.4 


160 


1.100 + 0.255 


4,15 


GX 17+2 


11.2+ 4.0 


750 


1.000 ± 0.250 


3.6 + 0.3 


650 


1.000 + 0.250 


4,16 


GX 340+0 


7.6+ 2.4 


500 


1.600 ± 0.400 


2.8 + 0.8 


350 


1.150 + 0.275 


4,17 


Cyg X-2 








3.6 + 0.6 


520 


0.800 ± 0.200 


4,18 


Sco X-l 


9.0+ 0.2 


780 


1.500 ± 0.375 


> 1.4 


590 


1.500 + 0.375 


4,19 



NOTES — The 1-<t errors are given. An aribtrary 25% error is used for the luminosity (see text) 
Qmax f° r 4U 0614+09 has been measured in the 4.6 — 60 keV energy range. 

The parameters for GX 340+0 and Cyg X-2 have been measured in the 5 — 60 keV energy range. 
The parameters for GX 17+2 have been measured in the 5.5 — 60 keV energy range. 
The parameters for Sco X-l have been measured over a variable energy range (see text). 
All the other parameter s have been measur e d in t he ful l PCA band. ^^^^^^^ 
REFERENCES — (1) Ivan Straaten et alJ teOOd): ( 21 Ivan Straaten et alj i2002h: (3 1 iFord et alJ <1997l): (4) 
iFord et alj fcOOfJl: (51 Thi s pape r; (6) IBarret et aI.U2005al): (71 lMendez et alj fcOOlft: ( 8) Ivan Straaten et al 
|2003|) : (91 bi Salvo et"aT| J2003ll : (10 1 iBarret et alj fe005bt) : (111 iBarret et alj |2005cJ)_; (121 |Pi Salvo et al 
1200 



120031 : (10 1 IBarret et alj J2005bll : (111 IBarret et alj |2005i., . , 

... (131 IBarret et alj i2005dY: (141 IFord et alj Il99^ : ( 151 Ijonker et alj 1200211: (161 iHoman et alj <20olk 
(171 Ijonkeret^rt . U2000I1 : (18l lWiinands et alj ^Uj98aTi ; (19) Ivan der Klis et alj Jl997t) 



of a source are not available in the literature; in some of 
those cases I myself have systematically analysed an pub- 
lished other QPO parameters (e.g., QPO frequency) of those 
sources, and since I still have the data available, I analyse 
them myself for this paper to obtain those missing QPO 
parameters. I describe those cases individually, including a 
summary of the analysis procedures, below. 

In Tables and [5] I provide all the measurements of the 
maximum coherence, Q max , and maximum rms amplitude, 
r max , for each QPO for each source, including the frequency 
of the corresponding QPO at which Q max and r max occur; 
in those Tables I also provide references to the papers from 
which these values were collected. 



2.2 Measurements over different energy bands 

In most cases in the literature, the QPO properties reported 
are those measured over the full energy band of the PCA; 
nominally, this range is 2 — 60 keV. In any case, the PCA 
instrument is not very sensitive above ~ 30 keV because 
of the lower effective area at those energies, and because 
the background dominates the signal. In a few cases (see 
Tables and |5J measurements were done over a different 
energy band to increase the sensitivity to kHz QPOs; since 
the rms spectrum of the QP Os increases steeply with energy 
(see e.g.. lBerger et al.ll996l1 . in those other cases the authors 
analysed the high-energy part of the data, above ~ 4—5 keV. 
On the other hand, the ageing of the PCA units (PCU), and 
changes applied to the gain voltage over the years to preserve 
the detectors affected the energy scale of the PCA as well as 
its effective area as a function of energy. The largest impact 
occurred at the low-energy end, since the low-energy cut-off 
of the detectors is a strong function of the instrument's high- 



voltage; from 1996 until present the lower energy boundary 
of the PCA increased by ~ 70% due to the ageing of the 
instrument and the gain changes applied to it 1 . 

Because of the strong dependence of the rms spectrum 
of the QPOs on energy, it is not straightforward to compare 
r measurements over different energy bands: The higher the 
low-energy boundary of the band over which one measures 
r, the higher the value of r. On the other hand, Q measures 
the number of cycles of the oscillation which, at least to first 
order, does not depend on energy. Therefore, one expects 
that the values of Q would not be much affected by the 
choice of the energy band over which they are measured. 

To confirm the above, and to assess the effect of 
the change of the energy band on the parameters, I pro- 
ceed as follows: I use the data for 4U 0614+09 from 
Ivan Straaten et alJ i2002i) who measured Q and r of both 
kH z QPOs over the full PCA band (2 - 60 keV, nominally), 
and lvan Straaten et alJ a200Cl) who measured the same prop- 
erties in the 4.6 — 60 keV range. In both cases they report 
Q and r as a function of QPO frequency, therefore I com- 
pare Q and r in both papers at the same QPO frequency. I 
find that Q is consistent with being the same in both energy 
bands. On the other hand, the rms fractional amplitudes in 
the 4.6 — 60 keV band are 1.12 + 0.04 (1-a error) times higher 
than the rms fractional amplitudes in the full PCA band. 

However, the energy spectrum of 4U 0614+09 is 
rather hard, and it may be that in other sources that 
have softer spectra, the rms amplitudes in different 
energy bands behave differently. To check this, I com- 
piled rms amplitudes in the 2 — 60 band for 5 sources: 



1 see http://heasarc.gsfc.nasa.gov/docs/xte/e-c_table.html 
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Table 2. Maximum coherence and maximum rms amplitude of the upper kHz QPO for the sources studied in 
this paper. 



Source 



(Hz) 



(L/L Edd )Qa 



(%) 



v n (L/L Edd ) r 
(Hz) 



Rcf 



4U 0614+09 


10.9 + 


2.8 


1270 


0.011 + 0.003 


16.5 + 0.3 


510 


0.005 ± 0.001 


1,2,3,4 


4U 1608-52 


12.5 + 


2.8 


1060 


0.050 + 0.012 


18.5 + 0.7 


475 


0.080 ± 0.020 


4,5,6,7,8,9 


4U 1636-53 


52.0 + 


14.0 


1230 


0.120 + 0.030 


17.1 + 0.6 


530 


0.028 ± 0.007 


4,5,10,11,12 


4U 1728-34 


10.5 + 


1.8 


1140 


0.070 + 0.017 


13.3 + 0.7 


600 


0.027 + 0.007 


2,4,6,13 


4U 1735-44 


















4U 1820-30 


23.0 + 


5.0 


1050 


0.250 + 0.062 


10.1 + 0.3 


580 


0.125 ± 0.031 


4,5,14,15 


Aql X-l 


















GX 5-1 


6.9 + 


0.8 


850 


1.600 ± 0.400 


2.9 + 0.2 


560 


1.100 + 0.275 


4,16 


GX 17+2 


13.0 + 


3.5 


1010 


1.000 + 0.250 


5.8 + 0.4 


590 


0.900 + 0.225 


4,17 


GX 340+0 


9.7 + 


2.8 


850 


1.600 + 0.400 


4.2 + 1.0 


530 


1.100 + 0.275 


4,18 


Cyg X-2 


22.2 + 


13.6 


960 


0.800 ± 0.200 


3.5 + 0.4 


860 


0.800 + 0.200 


4,19,20 


Sco X-l 


15.3 + 


0.3 


1040 


1.500 ± 0.375 


> 2.4 


875 


1.500 + 0.375 


4,21 



NOTES — The 1-<t errors are given. An aribtrary 25% error is used for the luminosity (see text) 
The parameters for GX 340+0 and Cyg X-2 have been measured in the 5 — 60 keV energy range. 
The parameters for GX 17+2 have been measured in the 5.5 — 60 keV energy range. 
The parameters for Sco X-l have been measured over a variable energy range (see text). 

All the other parameter s have been measur e d in t he ful l PCA band. 

REFERENCES — (1) Ivan Straaten et alj |20QO|>: (21 Ivan Straaten et alj fcOOSft: (3 1 iFord et alj Il997t>: (41 
iFord et al.l fcOOOft: (51 This pape r; (61 iMendez et al.l feOOlU; (71 iBarret et alj J2005al): (81 Ivan Straaten et all 

T^3T "\ In:..,".,.: 0_ T~» J l™r,<T" /1 r,\ IH^T^^TTTI F^^T. /TV _l I E^^T" (ml 



<2003l): 



(91 iGierlirlski fc Done! l2002f>: (1 1 |Pi Salvo et al 
lAltamirano et alJ <200fift : (131 iDi Salvo et alj ]200l|) : (141 " 
(161 iJonker et alj l2002t> : (1 71 iHoman et alj i2002ti: (181 
IWiinands et alj Jl998aF T(211 Ivan der Klis et alj Il997h 



20031): (111 IBarret et alj <2005t ): 
Altamirano et alJ 120051) : (151 iBloser et al, l200nll : 
XnkeretalJ j200Cl) : (191 lKuznetsoTT2002l) ; (201 



find that the rms fractional 
very well correlated (the 



4U 1608-52 jMendez et al] l200ll: Ivan Straaten et all 
20051). 4U 172 8 -34 jMendez fc van der Klis! Il999t 
" van Straaten et alJ 120021). 4U 1636-53 iBarret et all 
2005bJ ), 4U 0614+09 llvan Straaten et alJ EoOOl). an d 411 
1820-30 ijAltamirano et al J 1200.4 iBelloni et al.ll2005ll . and 
the co rresponding rms a mplitudes in the 5 — 60 keV 
range Jjonker et aljl200lfl . I 
amplitude in both bands are 
correlation coefficient is 0.97), and that the rms amplitude 
in the 5 — 60 keV band in these sources is 1.25 ± 0.02 times 
larger than the rms amplitude in the 2 — 60 keV. This 
number is slightly higher, but still consistent (within 3-a) 
with the value obtained above for 4U 0614+09. 



In summary, the above results show that: (i) One can 
compare directly Q-values in different energy bands, at least 
as long as the bands are not too different from each other. 
Because of this, in the rest of the paper I will compare Q 
values even if they come from different bands; (ii) changes 
of the low-end of the energy band affect the rms amplitudes 
by ~ 10 — 25% of its value or less (in the case of instrument 
gain changes and ageing, since changes of the low-energy 
boundary are smaller than those in the case of 4U 0614+09 
previously discussed, the effect is probably much less impor- 
tant). In the next section, when I plot the QPO parameters 
against luminosity and against each other, I divide the rms 
amplitudes in the band going from ~ 5 to 60 keV by 1.25 to 
make them compatible with the rms amplitude values ob- 
tained in the full PCA band. The values quoted in Tables 
and |21 however, are the ones actually measured and pre- 
sented in the original papers, in the energy band indicated 
there. 



2.3 Luminosity measurements 

Once I have collected all Q ma x and r max values for each QPO 
for each source, and the frequencies vg an d v u at which those 
maximum values occur, I use Figure 1 in lFord et al] (|200(J) 
to get the corresponding source luminosity: Using the QPO 
frequency as input, I read off the luminosity of the source at 
that frequency from that Figure. For three sources, 4U 1608- 
52, 4U 1636-53, and 4U 1820-30 (see below), the maximum 
rms amplitudes of L u occur in states of the source in which 
the QPO fr equencies a r e outs ide the range of frequencies in 
Figure 1 of Fo rd et al] l|2000j). In those three cases I either 
search the literature for flux measurements of the source 
in that state, or I extract spectra from the corresponding 
RXTE observations, and calculate the luminosities myself, 
or both. Below I explain those three cases in more detail. 

The luminosities that I report here are uncertain for a 
number of reasons: First, there is a statistical error in the 
fluxes derived from model fitting. Given that these sources 
are quite bright when kHz QPOs are detected, and that 
RXTE has a large effective area, these errors are small, usu- 
ally less than 5 — 10%. The second source of error is the 
accuracy with wh ich one can determ ine the luminosity us- 
ing the results of IFord et al] <|2000T) for a given QPO fre- 
quency. This is not so much how accurately one can read 
L from their plot, but the fact that for each source there 
is not a single X-ray intensity or X-ray flux value that cor- 
responds to a given QPO freq uency. This is the s o-called 
parallel-track phenomenon (e.g.. lMendez et alll999T) . In the 
extreme cases, the X-ray intensity of the source at a given 
QPO frequency may by as much as a factor of 3 different 
(e.g. 4U 1608-52. lMe"ndez et alll999h . On average, however, 
the range of intensitie s at fixed QPO frequency is smaller, of 
the order of 20-50% jMendez et all200ll) . The third source 
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of error is the use of the 2 — 50 keV flux as a measure of the 
bolometric flux o f the source. This is discussed in detail by 
iFord et alJ l)2000F) . and I therefore refer the reader to that 
paper. This effect contributes uncertainties of the order of 
20—25% of the reported luminosity values. The fourth source 
of error, if Q max and r max depended not on the luminosity 
of the source, L, but on the luminosity of the source normal- 
ized to its own Eddington luminosity, L/L-^dd, would be the 
use of a single Eddington luminosity to normalise the ob- 
se rved lumino s ities o f all the sources in the sample, as done 
bv lFord et alJ fgOOO) . Of course, if Q max and r max depended 
on L, and no t on L/J/Edd, this would be of no concern. The 
value used bv lFord et all ll200fj) is L E dd = 2.5 x f0 38 erg s _1 , 
which is the Eddington luminosity for a 1.9Mq neutron star 
accreting matter with cosmic abundance. Since all these sys- 
tems are thought to have had a common evolutionary path, 
and to have accreted between 0.4 and 0.8 Mq, the use of a 
single Eddington luminosity introduces an uncertainty that 
is of the order of 25% or less. (Notice that changes of the 
chemical composition of the accreting material may intro- 
duce a difference of up to a factor of 1.7 in the value of 
the Eddington luminosity; in the unlikely case that the neu- 
tron stars in this sample had different equations of state, 
the different redshifts at the surface of the neutron star may 
add an extra factor 1.5.) Finally, the largest uncertainty in 
the luminosity comes from t he error in the es timate of the 
distance to these sources fsee lFord et al.ll200ft for a discus- 
sion in the context of these sources); this uncertainty can 
introduce errors of up to 60% in the luminosity (see, e.g., 
IChristian fc SwankHl99?]). It should be mentione d that for 
GX340+0 and GX 5-l. lChristian fc Swankl lll997l) only give 
upper limits to th eir distance s , henc e for those two sources 
the luminosities in lFord et, alJ l)2000F) should be taken as up- 
per limits. 

In Tables [T] and |U and in Figure 01 use a fixed error of 
25% in the values of L/L^dd as indicative of the error in the 
luminosity. It is clear from the previous discussion that this 
is a lower limit to the real error in this quantity. 

2.4 Individual sources 

I now add a few remarks about special situations of some of 
the data for sources I include in this paper: 

4U 1608—52 : To m easure Q„ lax , I use the same data 
in lMendez et alJ (|200ll) . I accumulate power density spectra 
every 64 seconds of data up to a Nyquist frequency of 2048 
Hz; I average a variable number of consecutive power spec- 
tra, in all cases less than 15, in order to obtain significant 
QPO detections, and I search for high-frequency (> 250 Hz) 
peaks. If I detect two QPO peaks, the one with the lowest 
frequency is Le\ if I only detect one peak, I decide whether 
it is Li or L u on the basi s of the the rms-frequency relation 
(e.g., iMendez et alll200lfl . If it is L u I discard it, if it is Lg 
I measure Qe- 

Individual Qe measurements have relatively large errors. 
Since Qi is a function vg, I then sort the results on the basis 
of ve, and I calculate average Qe values over a narrow range 
of vi. As usual, this procedure reduces the errors in Qe at 
the expense of slightly undere stimating its maximu m value. 

I do not use the results of iBarret et alJ (l2005af) because 
in that case they calculate Q using whole observations; al- 



though th ey shift the individua l power spectra of each ob- 
servation iMendez et al.lfl998aL see below), still in many of 
those observations the frequency of the QPO drifts by sev- 
eral tens of Hz, which has a negative effect in the determi- 
nation of Q. Nevertheless, I find a Qf nax value that is larger 
than, b ut still consistent with, the value in IBarret et alJ 
J2005al) . 

For this source, rj,„ comes from Ivan Straaten et alJ 
(2003); they find this rms amplitude at a QPO frequency 
v u (see Table HI that is outside the frequency range in 
IFord et alJ J2000f) . and he nce no simultane ous luminosity 
measurement is available in lFord et alJlj200ffl . Van Straaten 
(2003) measured this rms amplitude duri ng their interval C; 
this c orresponds to regions 1 and 2 in iGierliriski fc Dong 
(2002), who find that the bolometric flux of the source in 
those regions is fbo i = 1 x 10 _9 erg cm -2 s _1 . Interval C of 
Ivan Straaten et alJ ll2003l) corresponds to the RXTE Obsld 
30062-01-01-05. I therefore downloaded those data from the 
RXTE archive, I produce a light curve of this observation 
and find that there are no X-ray bursts or any anomaly in 
the data; I then extract an X-ray spectrum using the Stan- 
dard2 data following the procedures described on the RXTE 
web pages 2 . I fit this spectrum with a model consisting of a 
blackbody and a power law plus a Gaussian emission line at 
around 6.5 keV, all affected by interstellar absorption. The 
model fits the data well, with a \ 2 P er degree of freedom 
of about 1. The 2 — 50 keV unabsorbed flux of the source 
in this obse rvation is 1 x 10~ 9 erg cm " 2 s _1 , consistent with 
the flux of iGierlihski fc Donej j|2002j). Using the same dis- 
tance to this source as in TFordet^ll ll200fl ). d — 3.6 Kpc, 
this flux corresponds to a luminosity L — 1.96 x 10 37 erg s _1 , 
or L/Z/Edd = 0.0 8 using the Eddington luminosity used by 
IFord et alJ l)2000l) . L E dd = 2.5 x 10 38 erg s" 1 . 

4U 163 6—53: To measure Q max; I use t ne same data pre- 
sented in iDi Salvo et all (120031) . The procedure to produce 
power spectra, identify the QPOs, and measure Qe is the 
same as the one described for 4U 1608-52. The value of 
Q max that I find is slig htly higher than, but still consistent 
with, the one found by IBarret et~ai] ll2005bl). 

For this source, r^aa comes from lAltamirano et alJ 
( 2006); they find this rms amplitude at a QPO frequency u u 
(see Tabled that, as in the case of 4U 1608-52, is outside the 
frequency range in lFord et alJ (I2000T) . and henc e no simul- 
taneou s luminosity measu r ement is available in IFord et"al] 
l|2000l) . lAltamirano et al"l (l200rJ) find this rms amplitude 
value during their interval C, which corresponds to the 
RXTE Obsld 60032-05-10-000 and 90409-01-01-02. I pro- 
duce a light curve of these observations and find that there 
is an X-ray burst in the first one, and no X-ray bursts or 
any anomaly in the data of the second one; I then follow the 
procedures described on the RXTE web pages to extract an 
X-ray spectrum from the Standard2 data of the second ob- 
servation. I fit this spectrum with a model consisting of a 
blackbody and a power law plus a Gaussian emission line at 
around 6.5 keV, all affected by interstellar absorption. The 
model fits the data well, with a \ 2 P er degree of freedom of 
about 1. The 2 — 50 keV unabsorbed flux of the source in 
this observation is 1.9 x 10 ~ 9 erg cm" 2 s -1 . U sing the same 
distance to this source as in lFord et alJ J200CT) . d = 5.5 Kpc, 



2 http:/ /heasarc. gsfc.nasa.gov/docs/xte 



Amplitude and coherence of the kHz QPOs 7 




Figure 3. Left panel: The maximum coherence (upper panel) and maximum rms amplitude (lower panel) of the lower kHz QPO for the 
sources listed in TableQ] as a function of the source luminosity at the time at which those maximum values were reached. The luminosity 
is in units of the Eddington luminosity for a 1.9Mq neutron star. Filled circles indicate measurements over the full energy band covered 
by the PCA on board RXTE. Open circles indicate measurements over a limited energy band (see Table 111 for details); in these cases 
(except for Sco X-l), the rms amplitudes have been divided by 1.25 (see text for details). The rms amplitude in the case of Sco X-l 
is not corrected for dead-time, and hence is only a lower limit (indicated with an arrow pointing upwards in the lower panel). Right 
panel: The maximum coherence (upper panel) and maximum rms amplitude (lower panel) of the upper kHz QPO for the sources listed 
in Table |21 as a function of the source luminosity at the time at which those maximum values were reached. The luminosity is in units 
of the Eddington luminosity for a 1.9Mq neutron star. Symbols are the same as in the left panel. 



this flux corresponds to a luminosity L = 6.9 x 10 36 erg s 1 , 
or L/L^dd = 0.0 28 using the Eddington luminosity used by 



iFord et alJ d200d) . L Ed d = 2.5 x 10 38 erg s" 1 



and Q" 



I use the 



4U 1820-30: To mea s ure ( 
same data as in [Mendej d2002 



and iBelloni et alJ (2005). 
The procedure to produce power spectra, isolate Le, and 
measure Qe is the same as the one described for 4U 1608-52 
and 4U 1636-53. To detect L u and measure Q u , I proceed 
as follows: I group the power spectra on the basis of the fre- 
quency of Le, ve, such that within each group ve does not 
change by more than 50 Hz. I then shift the power spec- 
tra such that in each group the frequency of Le is the same 
in all power spectra, and then I calculate an average power 
spectrum per group. This procedure eliminates the effect 
of the drift in Le, and since the frequency diff erence in this 
source is more or less constant when vt changes dZhang et alJ 
Il998al) . it also corrects the effect of the drift in L u . I then fit 
the power spectra with two Lorentzians that represent the 
QPOs, and a constant that represents the Poisson counting 
noise. The fit yields, among other parame ters, Q u . 

For this source, the rjjj a:c value is from lAltamirano et alJ 
(2005); they find this rms amplitude at a QPO frequency 
v n (see Table |5J that, as in the case of 4U 1 608-52 and 
4U 16 36-53, is outside the frequency range in IFord et ail 
(2000), and he nce no simu l taneou s luminosity measu r ement 
is available in IFord et alJ d2000l) . lAltamirano et alJ (120051) 
report this rms amplitude at a position in the colour-colour 
diagram of 4U 1820-30 that is consistent with the interval 
1 of lBloser et ail (|2000). {S a is a parameter that mea- 
sures the position along the C-shaped track traced out by 
the sou rce in a colour-col o ur dia gram.) During their obser- 
vation, lAltamirano et all (120051) find that the intensity of 
the s ource is 297 counts s" 1 PCU -1 , while in their observa- 
tion, [Blose^^^ <l2000fi find the source at 312 c ounts s _1 
PCU" 1 . For their observation. iBloser et ail d2000l) find a 2- 
50 keV luminosity of 2.25 — 2.31 x 10 37 erg s _1 for an assumed 



distance d — 6.4 Kpc. IFord et"al] <l2000l) use a distance to 
NGC 6624, the globular cluster that contains 4U 1820-30, 
d = 7.5 Kpc. To be consistent with the other luminositie s 
for this source, which were taken from IFord et alJ J2000^ . 



I convert the luminosity in IBloser et al 



2000) to the one 



corres ponding to the distance to 4U 1820-30 in lFord et alJ 
(2000). This corresponds to L = 3.15 x 10 37 erg s~\ or 
L/L^dd = 0.125 using the Eddington luminosity used by 
IFord et alJ fcOOd) . L Edd = 2.5 x 10 38 erg s -1 . 

4U 172 8—34: To mea s ure Q may , I use the same d a ta pre - 
sented in |Pi Salvo et alJ d200lj) and IMendez et al.l fcOOll) . 
The procedure to produce power spectra, identify the QPOs, 
and measure Qe is the same as the one described for 4U 
1608-52. 

4U 1735-44: For this source, iBarret et all d2005cl) re- 
cently reported a systematic study of Qe and re as a function 
of vt ; I therefore are reported there. 

There are, howev er, very few reports of the properties of L u 
in the literature dWiinands et alJll998bt IFord et alJll998t) : 
no reliable value of Q ma x or r max are available (see also Fig. 
. I therefore do not discuss L u of this source in this paper. 

Aql X-l: For this source only one kHz Q PO has been de- 
tected dZhang et al.lll998bl:ICui et alJll998l) . which appears 
to be Le IMendez et alJl20 01*l. and hence there are no QPO 
parameters available for L u . 

Cyg X-2: Since there are only three measurements of Qe 
for Cyg X-2 (Fig.0, I do not consider Qe for this source in 
the rest of this paper. Also, there are only t wo measurements 
of re f or this source, but the upper limits dWiinands et alJ 
1998a) seem to indicate that the maximum r^-value has been 
measured, and hence here I use re in the analysis. 

Sco X-l: To avoid detector safety triggers, telemetry 
saturation, and to reduce the dead-time effects produced 
by the high count rate of Sco X-l, some observations of 
this source were carried out with the source slightly off-axis, 
with some of the five proportional counter units of the PCA 
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Figure 4. The maximum coherence (upper panels) and maximum rms amplitude (lower panels) of the lower kHz QPO (left) and 
the upper kHz QPO (right) for the sources listed in Tables and [5] as a function of the average sou r ce hardness, HI, defined as 
the r atio of the 40 — 80 keV to the 13 — 25 keV count rate, measured with HEAO-1 jLevine et al, 1984; van Paradii s fc van der Klisl 
1994). For an explanation of the symbols see the caption of Figure[3] The rms amplitudes in the case of Sco X-l are not corrected for 
dead-time, and hence are lower limits (indicated with a vertical arrow in the lower panels). Upper limits to the hardness are indicated 
with horizontal arrows pointing to the left. I do not include GX 5-1 in this Figure because the HEAO-1 measurem ents of this source 
suffer cont amination from a pr eviously unknown hard X-ray source in the field (see explanation in lLevine et al.l|l984fl . most likely GRS 
1758-258 iGilfanov et alJll993D . 



switched off, recording only photons detected by the up- 
per anode chain of the PCA, recording only photons from a 
limited e nergy range, or using a c o mbination of these con- 
strain ts II van der Klis et al.l Il997t IMendez fc van der Klisl 
2000). Despite all these efforts, the source count rate during 
these observations remained high, at a level in which the 
dead-time of the PCA is still unknown. This m eans that the 
rms values (taken from lvan der Klis et al.ll997t) are very un- 
certain, and in fact could be larger than reported. For that 
reason, for Sco X-l in this paper I report the maximum rms 
amplitude as a lower limit. 



3 RESULTS: 

Figure |H] shows the dependence of the maximum coher- 
ence and maximum rms amplitude of both kHz QPOs as 
a function of source luminosity in units of the Edding- 
ton luminosity for a 1.9M0 neutron star. From this Fig- 
ure it is apparent that r^ ax and r^ ax both decrease more 
or less exponentially with L/L-Edd- For rf nax the e- folding 
scale is L/L-^dd = 0.86 ± 0.06, while rj^ ax decreases signifi- 
cantly faster with L/Lsdd than r max , with an e- folding scale 
L/LiEdd = 0-61 ± 0.0 3 . Usin g roughly the same sample of 
sources, Ijonker et al.l (1200 ll) had already noticed that the 
rms amplitude of the upper kHz QPO decreases as the lu- 
minosity of the source increases. 

From the same Figure, it is also apparent that at low lu- 
minosity, Qmax nrs t increases with L up to L/Lsdd ~ 0.04, 
and then it decreases exponentially, at a faster rate than 
r max> with an e-folding scale L/Z/Edd = 0.47 ± 0.01. On the 
other hand, Q^ax does not show any significant trend with 
luminosity; although there is a hint that it stays constant up 
to L/L-Edd ~ 0.07, then it increases up to L/Z/Edd ~ 0.12, 
and after that it decreases for increasing luminosity, statis- 
tically a fit with a function that represents that behavior is 
not significantly better than a fit of just a constant. 



The gap in Figure EH at L/Z/Edd ~ 0.25 - 0.7 sep- 
arates the Atol l sources at low L, and the Z sources 
at high L (see lHasinger fc van der Klisl [l989, for a def- 
inition of Atoll and Z sources). That gap would be 
occupied by the four intermediate-type sources, GX 
9+1, GX 9+9, GX 3+1, and GX 13+1, which so 
far have not shown any kHz QPOs (ISjro hjnayer] [l99g|; 
Wiinaiidsjar^terlCh^L^aii^aradiis 1998 a iHoman et alJ 
ll99SUSchnerr et alJl2003l : lOosterbroek et al.ll200ll) . The ~ 
per limit to the rms amplitude of the QPO in these sources 
ranges from 1.6% to 2.6%. S ince the range of lumino si- 
ties spanned by these sources iChristian fc Swanklll997l) is 
L/L-Edd ~ 0.12 — 0.44, these upper limits are much lower 
than would be expected from the interpolation of the trends 
of rJn ax and r^ ax with L/L E dd in Figure [3 

From Figure [3] it is apparent that the maximum rms 
amplitude of both kHz QPOs and the maximum coher- 
ence of the lower kHz QPO are consistently lower in the 
Z sources than in the Atoll sources. A valid question would 
then be whether this difference could be due to a bias in 
the way in which the QPOs of Z and Atoll sources are mea- 
sured, combined with the frequency drift of the QPOs (see 
H2.H . A comparison, however, of the width of the lower kHz 
QPO over a short time interval in a Z and an Atoll source 
shows that this is not the case (as I describe in §2, the 
rms amplitudes are less affected by the frequency drift). To 
see this, I compare an RXTE ob servation of the Z source 
Sco X-l from May 25 1996 (see Ivan der Klis et alJ Il997l : 
IMendez fc van der Klisl lioOO). and another RXTE observa- 
tion of the atoll sou rce 4U 1608-52 from March 26 1998 (see 
IMendez et ai] [l999). In the Sco X-l observation the power 
spectrum shows two simultaneous kHz QPOs that move in 
frequency; Li moves from vi ~ 590 Hz to vg ~ 650 Hz, 
and L u moves from v u ~ 890 Hz to v u ~ 950 Hz. In the 
observation of 4U 1608-52 the power spectrum also shows 
two simultaneous kHz QPOs that move in frequency. In this 
case Lg_ moves from vi ~ 590 Hz to vt ~ 610 Hz, and L u 
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Figure 5. The maximum coherence of the kHz QPOs as a func- 
tion of the maximum rms amplitude of the same kHz QPO, for 
the sources listed in Tables and |5] the upper and lower panel 
show the upper and the lower kHz QPO, respectively. For an ex- 
planation of the symbols see the caption of Figure[3] In the case of 
Sco X-l the rms amplitudes are not corrected for dead-time, and 
hence are lower limits (indicated with horizontal arrows pointing 
to the right). 



moves from v u ~ 890 Hz to i/ u ~ 910 Hz. To compare the 
QPOs in both observations over the same frequency range, 
I calculate an average power spectrum of the Sco X-l data 
using only intervals in which 590 Hz ^ vt ^ 610 Hz; this 
corresponds to 24 individual power spectra each of them 
16-s long. In the case of 4U 1608-52 I calculate an average 
power spectrum of the whole observation, corresponding to 
41 individual power spectra each of them 64-s long. To cor- 
rect for any residual drift in the QP O frequency, I apply th e 
shift-and-add procedure described in lMendez et alJ dl998aTl . 
The results, however, do not change significantly if I aver- 
age the individual power spectra without first applying this 
method, most likely because the frequency drift is small in 
both cases. From the above power spectra I find that for Sco 
X-l Q e = 4.2 ±0.4, whereas for 4U 1608-52 Q e = 74.0 ±4.6, 
which shows that in the Z source Sco X-l, even over very 
short time intervals (in this case 384 s), the lower kHz QPO 
is significantly broader than in the Atoll source 4U 1608- 
52. This in turn shows that the drop of the maximum QPO 
coherence of the lower kHz QPOs at high luminosities in 
Figure |3] is real. 

Tabled shows that in all sources the maximum coher- 
ence factor of the lower kHz QPO, Qmax, occurs more or 
less at the same frequency (the spread is ~ 17%); Table |2] 
shows that the same is true for the maximum coherence fac- 
tor of the upper kHz QPO, Qn lax (the spread in this case 
is ~ 13%). Since the lifetime of an oscillation of frequency 
v and coherence factor Q is roughly r ~ Q/u, the results 
described above also indicate that the maximum lifetime of 
Lt, T"max! nrs t increases and then decreases with L/L^dd, 
whereas the maximum lifetime of L u , T max , is independent 
of L/Z/Edd- 



IVan Paradiis fc van der Klij dl994f) have shown that 
there is a general correlation between the average source lu- 
minosity and the average source hardness, HI, defined as the 
ratio of the count rate in the 40 — 80 keV band to the count 
rate in the 13 — 25 keV band, measured with HEAO-1 (see 
iLevine et al]ll984h . In that respect, as expected, the plots 
of Qmax and r max vs. hardness in Figure 2] show that r max 
and r^ ax both increase with the hardness ratio HI, Qm ax 
increases with HI and then it decreases for 4U 0614±09, the 
hardest source in this sample, and as in the plots as a func- 
tion of luminosity, Qm ax is consistent with being constant 
with HI. At the risk of pointing out something obvious, the 
opposite behaviour of the coherence and amplitude of the 
QPOs to that in Figure [3] is due to the anticorrelation be- 
tween HI and L/L^dd- I caution the reader that contrary 
to the L/L-Edd values that I present in Tables and [5] and 
I plot in Figure [31 for this Figure I use average values of 
the spectral hardness measured several years before the kHz 
QPOs were discovered. Notice that GX 340±0, for which the 
distance, and hence the luminosity, is uncertain (see § 12.31 . 
appears in Figure 0] close to the other Z sources (open sym- 
bols; GX 340±0 is the point at HI = 0.36, just to the left of 
GX 17+2 at HI = 0.37); since HI is a distance-independent 
parameter, this suggests that the distance to GX 340+0 is 
not too much in error. 

It is interesting to notice in this Figure that there is 
no gap in the distribution of sources as a function of hard- 
ness between the Z and Atoll sources. This is opposite to 
what is apparent in the plot of QPO parameters vs. luminos- 
ity, where there is a gap corresponding to the intermediate- 
type, the GX-sources (see above). The HI values of the GX- 
sources range from 0.03 to 0.14. However, both in Figure 
and 2] there appears to be a gap in the values of Q^ ax , and 
perhaps also r max and rj,,, between the two types of sources 
(see also below). 

Figure [5] shows the relation between Q max and r max for 
the same QPO. To produce this Figure (and the next one), 
I plot the Q max vs. the r max values for the same source, 
even if they occur at slightly different luminosities. I also 
produced plots of Q max vs. r max combining values for sources 
that are closest in luminosity, even if this means combining 
values from different sources, but the trends described below 
remain the same as with the previous choice. 

Except for the case of 4U 0614+09, Qm ax an d J"m ax 
are positively correlated with each other. Concerning 4U 
0614+09, the hardest source in the sample and the one at the 
lowest luminosity, it is as if in this case Q^ax were too low for 
r m aX i or conversely, as if r^ax were too high for Q^ ax . This 
is already apparent in Figure [3] where Lt of 4U 0614+09 
(low luminosity part of the plots) looses coherence without 
a similar decrease of its rms amplitude. Figure [5] also shows 
that Qm ax is independent of r^ ax . 

Let us examine the case of the low Q^ ax value in 4U 
0614+09 to check whether it is real: The measureme nts of 
Qi in this source come from Ivan Straaten et alJ i2000t) . and 
are made over a different energ y band than those of the 
other sources. For 4U 0614+09 Ivan Straaten et alJ feOOfJ) 
use a band that starts at 4.6 keV, whereas for other sources 
at slightly higher luminosity the Qe values are measured 
using the whole PCA band, nominally starting at 2 keV (see 
Table0 . In H2.2I I compare measurements of Qe as a function 
of vi (and also Q u as a function of v u ) for this source over 
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Figure 6. The maximum coherence and maximum rms amplitude of one kHz QPO as a function of the maximum coherence and 
maximum rms amplitude of the other kHz QPO, for the sources listed in Tables □ and El The upper and lower left panels show plots, 
respectively, of Q m ax and r max vs. Q m ax! * ne upper and lower right panels show plots, respectively, of Q ma x and r max vs. r max . For an 
explanation of the symbols see the caption of Figure^ In the case of Sco X-l the rms amplitudes are not corrected for dead-time, and 
hence are lower limits (indicated with arrows pointing up and to the right). 



the full PCA band (2 - 60 keV: Ivan Straaten et ai1l2002fl 
and over the 4.6 — 60 keV band llvan^toaten^^dl 2000), 
and I find that both sets of values are compatible with each 
other. This indicates that most likely the difference in Q ma x 
between 4U 0614+09 and other sources at slightly higher 
luminosity is real. Recently, Barret (private communication) 
did a similar analysis on the data of 4U 0614+09 as he did 
for instance for 4U 1636-53 ijBarret et alJ l2005bTl . and he 
finds that Qe in 4U 0614+09 is ~ 6 times smaller than in 
4U 1636-53, consistent with the results in Table This 
confirms that Q ma x i n 4U 0614+09 is in fact smaller than in 
the other sources at similar but slightly higher luminosity, 
and hence that the decrease of Q ma x a t low luminosity is 
real. 

As I already mentioned, there appears to be a gap be- 
tween the Q max and r max values of the Z and Atoll sources, 
with the coherence showing the largest gap (Figure 3, lower 
panel). 4U 0614+09 appears to be the only source to (par- 
tially) break this rule, since it has a Q max value that is 
intermediate between those of Z and Atoll sources. Simi- 
larly, the upper panel of Figure 3 shows a gap between the 
r max values, but not for the Q ma x values, of the Z and Atoll 
sources. While this can indicate a dependence on luminos- 
ity (see Fig. [3J or on spectral hardness (see Fig. [IJ, this 
could also point to a difference between Z and atoll sources. 
It is worth noting, however, that there is still a trend of 
Q max and both rms amplitudes within the atoll sources in 
Figure [3] Furthermore, there is a significant trend in the 



relations between Q max and r max within the atoll sources; 
e.g., the relation between Q max an d r max (lower panel of 
Figure |SJ is 8a different from a constant. All this suggests 
that the distinction between Z and Atoll sources cannot be 
the (only) explanation for this difference. 

The plots of Q m ax and r max of one QPO vs. the same 
parameters of the other QPO are shown in Figure HJ From 
this Figure there are three things apparent: (i) Q max is inde- 
pendent both of Q max and r max (upper panels); (ii) as in the 
case of the plot of r max vs. Q max , except for 4U 0614+09, 
r max is positively correlated with Q max . As in the case of 
Le , it appears as in 4U 0614+09, the source with the lowest 
luminosity in the sample, the amplitude of L u , r max , were 
too high for the coherence of Le, Q max , or alternatively, the 
coherence of Le, Q max , were too low for the amplitude of 
L u , r max (lower left panel). Finally, (iii) the amplitudes of 
both kHz QPOs are positively correlated with each other, 
including in this case the low-luminosity source 4U 0614+09 
(lower right panel). 

From the lower panels of this Figure, it is also apparent 
that there is a gap between the r max values of the Z and 
Atoll sources, similar to the ones described above 

To summarize the results from Figures |S] and HJ the 
maximum amplitude and coherence of Le, r max and Q 



i 

max 5 



and the maximum rms amplitude of L u , 



are all corre- 



lated with each other (in the case of Q m axj at least above 
L ~ 0.04LEdd), whereas the coherence of the upper kHz 
QPO, Qmax, is independent of all the other parameters. 
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Figure 7. The rms amplitude (upper panels) and coherence (lower panels) of the lower kHz QPO (left panels) and the upper kHz QPO 
(right panels) for 4U 1608—52 as a function of the hardness of the source. The hard ness is defined as the logarithm of the ratio of the 
count rate in the 9.7 — 16.0 keV band to the count rates in the 6.4 — 9.7 band f see iMendez et alj|l999h . The lines that connects the 
points show the evolution of QPO frequencies, which generally increases from right to left in these plots. 



4 DISCUSSION 



I study the maximum amplitude, r max , and maximum co- 
herence, Q m ax, of the kHz QPOs as a function of luminosity 
and hardness for a large sample of low-mass X-ray binaries. I 
show, for the first time, that the maximum coherence of the 
lower kHz QPO, Qm ax , first increases up to L ~ 0.04LEdd 
and then decreases with luminosity, whereas the maximum 
coherence of the upper kHz QPO, QJ lax , is independent of 
luminosity. I also find that the maximum rms amplitudes 
of both the lower and the upper kHz QPOs, r^ax and r^ax, 
respectively, decrease monotonically with luminosity and in- 
crease monotonically with the hardness of the source. (The 
dependence of r^, ay on lumin osity and hardness was first 
reported bv l.Tonker et aljEoOll) . 



From the above results it follows that for all sources, 
r max an d rjjjax are positively correlated with each other. 
Also, for all sources with L > 0.04Z/Edd, that is all sources 
in this paper except 4U 0614+09, the hardest source in the 
sample and the one at the lowest luminosity, Q^ax is posi- 
tively correlated both with rf nax and r„ lax . QJ^ax is indepen- 
dent of Qmax °r the maximum rms amplitude of the kHz 
QPOs. 



Since the frequencies at which the values occur are 

more or less the same in all sources, and the same is true for 
Qm a x 5 although the frequencies in this case are higher, the 
dependence of Qmax an( A Qm ax on luminosity also reflects 
the dependence of the maximum lifetime of the QPOs on 
luminosity (and spectral hardness). Therefore, for the lower 
kHz QPO, the maximum QPO lifetime first increases and 
then decreases with luminosity. For the upper kHz QPO, the 
maximum QPO lifetime is independent of luminosity. 



4.1 The relation between the ISCO and the drop 
of QPO coherence and rms amplitude 

In individual sources, both r> and Qe increase with vt 
and then drop rather abruptly at the high end of the ui 
range; r u also increases and then drops at high v u val- 



with j/„ (e.g.. |Di Salvo et al 
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2006). In the case of 4U 1636-53, 
erpret the sudden drop of the coher- 
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ence and rms amplitude of Le, together with the existence 
of a frequency above which Le is not detected, as evidence of 
the innermost stable circular orbit, ISCO, around the neu- 
tron star in this system. In general relativity, different from 
the Newtonian theory of gravitation, the effective potential 
as a function of radial distance to the central source has a 
maximum. A particle in a circular orbit at that radius would 
be in unstable equilibrium; if perturbed, the particle would 
fall onto the central object. No stable orbit around the cen- 
tral object is possible inside the radius of the ISCO, which 
in the Schwarzsc hild case (non-rotating central obj ect) is 
nsco = 6GM/c 2 (iBardeen. Press, fc Teukolskvlll972l) . 

From the results of individual sources (see references 
above) and those of the sample of sources that I present in 
this paper, it is apparent that the behaviour of the coher- 
ence and rms amplitude of the kHz QPOs as a function of 
the QPO frequency in individual sources is similar to the be- 
haviour of the maximum coherence and maximum rms am- 
plitude of the kHz QPOs as a function of luminosity in the 
sample of sources. Since in individual sources there is a gen- 
eral relation between QPO frequencies and source intensity, 
in the sense that at higher intensity the QPOs generally ap- 
pear at higher frequencies (but remember the parallel-track 
effect), this raises the question of whether the same mecha- 
nism may be behind both behaviours. 

At first sight, this may appear problematic because it is 
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known that the relation between QPO frequency and inten- 
sity in individual sources is more complex than the gen- 
eral frequency-luminosity trend described above, and de- 
tailed QPO frequency vs. intensity pl ots show the so-calle d 
parallel-track phenomenon (see e.g.. IMendez et ail 1 1999). 
However, the link between these two behaviours need not 
be luminosity, but could be the high-energy emission (or 
hardness) in these systems. On one hand, in the sample of 
sources the maximum coherence and maximum rms ampli- 
tude of the kHz QPO, except the maximum rms amplitude 
of L u , appear to correlate fairly well with spectral hard- 
ness (see Figure 0J, while on the other hand in individ- 
ual sources the QPO frequencies are well correlated wit h 
the spectral hardness of the source iMendez et aljll999l) . 
the index of the power l aw that fits the hi gh-energy part 
of the X-ray spectrum ]Kaaret et al J 119981) . or S, a pa- 
rameter that measures the position along the track traced 
out by the source in a colour-colour or colour-intensity di- 
agram (called S z a nd S a in the Z and Atoll sources, re- 
spectively; see e.g., iHertz et al.l Il992l Ijonker et aD |2000; 
Men dez fc van der Klislll999ir From this, it follows that in 
individual sources there should be a relation between QPO 
rms amplitude and coherence on one hand and spectral hard- 
ness on the other. To my knowledge, a plot like that has 
never been published, but can be easily constructed using 
published results. For instance, in the case of 4U 1608-52, 
this can be done by combini ng the plots of QPO frequency 
vs. hardness of Figure 3 in Me ndez et al] <ll999l) and the 
plots of QPO coherence and QPO rms ampli tude vs. QPO 
fr equency of Figure 2 i n lBarret et al] (l2005cT) and Figure 3 
in IMendez et al.| ]200 J) , respectively. In Figure|7|I show the 
rms amplitude (upper panel) and coherence (lower panel) of 
the lower kHz QPO (left panel) and the upper kHz QPO 
(right panel) in 4U 1608-52 as a fun ction of hardness. In 
this case the hardness is defined (see Me ndez et alj|l999l) 
as the logarithm of the ratio of the PCA count rate in the 
9.7 - 16.0 keV band to the PCA count rates in the 6.4 - 9.7 
band. (Notice that the hardness in this plot is denned over 
a different energy band than in the case shown in Figure 0]) 
The lines connecting the data points indicate the evolution 
of QPO frequency, which generally increases from right to 
left. This Figure shows the drop of the rms amplitude and 
coherence of the lower kHz QPO and of the rms amplitude 
of the upper kHz Q PO at the very high end of the QPO fre- 
quency range (c f |Di^5cdvoet al l l200ll . l2003t IMendez et al] 
l200ll; (Barret et al]l2005allblc]) 

The trend seen in Figure |7| for the coherence of the 
lower kHz QPO in 4U 1608-52 is rather similar to the be- 
haviour of the maximum coherence of the lower kHz QPO 
for the sample of sources in Figure^Kleft panel; notice again 
that the hardness is defined over a different energy band in 
each Figure): When 4U 1608-52 is hard, to the right of Fig- 
ure [7] the coherence is low, while for the hardest source in 
the sample, to the right of the left panel of Figure 0] the 
maximum coherence is also low. When 4U 1608-52 becomes 
softer, to the left of Figure |7| the coherence first increases, 
and then drops abruptly (see black points in Figure |7] lower 
panel), while in the sample of sources the maximum coher- 
ence behaves in a similar way. The behaviour of the rms 
amplitude in 4U 1608-52 and the maximum rms amplitude 
in the sample of sources is more or less similar, except that 
in the case of 4U 1608-52 when the source is hard the rms 



amplitude is low, increases as the source gets softer and then 
drops abruptly (see black points in Figure |7| upper panel), 
whereas in the sample of sources moving from hard to soft 
sources the maximum rms amplitude remains more or less 
constant and at the end drops. To summarize, the relevant 
conclusion from this comparison is this: Both the rms am- 
plitude and coherence of the lower kHz QPO in 4U 1608-52 
and the maximum rms amplitude and maximum coherence 
of the lower kHz QPO in the sample of sources drop abruptly 
when the hardness decreases. 

The comparison between individual sources and the 
sample of sources suggests that the same mechanism is re- 
sponsible for the drop of coherence and rms amplitude of the 
lower kHz QPO with QPO frequency in individual sources 
as well as for the drop of maximum QPO coherence and 
maximum QPO rms amplitude with luminosity in the sam- 
ple of sources. The comparison between 4U 1608-52 and 
the sample of sources in the previous paragraph suggests 
that most likely the mechanism is related to the high-energy 
emission in these systems. This does not necessarily mean 
that the fractional emission at high energies (represented by 
the hardness or X-ray colors) is the root mechanism that 
drives all QPO parameters (QPO frequency, coherence, and 
rms amplitude). For instance, one possibility (there could 
be many others) is that the (instantaneous) m ass accretion 
rate s ets the size of the inner radius of the disc llvan der Klij 
2001), which in turn determines the QPO frequency, as well 
as the relative contribution of the high-energy part of the 
spectrum to the total luminosity. If the efficiency of the 
modulation mechanism and the lifetime of the oscillations 
that produce the QPO depended upon the emission from 
the high-energy part of the source spectrum (see § 14.21 for 
a discussion of possible ways in which this could happen), 
observationally it would appear as if the coherence and rms 
amplitude of the QPO depended upon the QPO frequency, 
and hence upon the radius in the disc at which the QPO is 
produced. The sudden drop of the coherence and rms am- 
plitude of the QPO at some QPO frequency would then be 
associated to a dynamical peculiarity in the accretion disc, 
for instance the ISCO. Observing the same source repeat- 
edly would not allow to distinguish the above scenario from 
one in which QPO coherence and rms amplitude were actu- 
ally set by QPO frequency or the dynamics in the accretion 
disc. 

To distinguish between these two options, one would 
need to observe a sample of sources of kHz QPOs for which 
the mass-accretion rate, and hence the relative contribution 
of the high-energy part of the spectrum to the total emis- 
sion, was different from source to source. In that case, QPO 
coherence and rms amplitude would drop for sources accret- 
ing mass at higher rates, even if the frequency of the QPO 
was more or less the same from source to source. Since, as I 
have shown in this paper, this is the general behaviour ob- 
served in sources of kHz QPOs, it is reasonable to infer that 
a mechanism similar to the one I described in the previous 
paragraph is effective in setting the coherence and rms am- 
plitude of the kHz QPOs. If this is correct, this also implies 
that the drop of QPO coherence and rms amplitude as a 
function of QPO frequency in individual sources cannot be 
due to effects of the ISCO. 

Note also that in individual sources not just the rms 
amplitude of the kHz QPOs, but also the rms amplitude of 
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other lower-frequency QPOs decrease with increasing QPO 
frequency. For instance, in four Atoll sources, 4U 1728-34, 
4U 1608-52, 4U 0614+09, and 4U 1636-53, the rms ampli- 
tudes of the "bump", a QPO at ~ 0.1 — 30 Hz, the "hump", 
a QPO at ~ 1 - 40, and the hectohertz QPO at ~ 100 - 300 
Hz, all drop as the frequencies of the kHz QPOs increase, in a 
similar fashi on as the amplitude of t he upper and lower kHz 
QPOs (e.g.. lAltamirano et al.ll2006t and references therein; 
see also there a description of these other QPOs). This also 
argues against the interpretation of the ISCO as the cause 
of the drop of the rms of the kHz QPOs, and indicates that 
the amplitudes of all variability components are set by the 
same mechanism which, as I suggested, could be the same 
one that governs the high-energy spectral component. 

The idea that the rms amplitude and coherence of the 
QPOs depend on the high-energy emission, and not on QPO 
frequency, can in principle be tested; for instance, if one 
of the kHz QPO sources showed a sudden change in the 
high-energy part of the spectrum, but the change of QPO 
frequency was less sudden. This could for example occur 
if the component that sets the high-energy emission in the 
spectrum of these sources and the one that sets the QPO 
frequency had different responses to changes of the mass ac- 
cr etion rate. This sce nario is similar to the one proposed 
bv Ivan der Klisl (2001) to explain the parallel tracks in the 
frequency vs. intensity plots of these sources. If this was 
the case, one would observe a kHz QPO that had a coher- 
ence and rms amplitude that would not match those ex- 
pected on the basis of th e Q — v and r — v relations in, e.g., 
iBarret e ^ilj jj^S^Hbloj ). I w ould like to point out that Fig- 
ure 1 in Barret et al] £2006) appears to show an effect like 
this: Individual measurements of the coherence of the lower 
kHz QPO at the same QPO frequency differ significantly 
from one another. E.g., in 4U 1636-53 at Vf ~ 820 Hz, Qi 
ranges from Q e = 80 ± 20 to Q e = 180 ± 10. This difference, 
however, may still be due to uncorrecte d drifts of th e QPO 
frequency during the intervals in which IBarret et al] (I2006T) 
measured the coherence of the QPO, although according to 
their description they had taken this into account in pro- 
ducing their Figure. 

4.2 The modulation mechanism and the lifetime 
of the kHz QPOs 

How can the mechanism that produces the high-energy emis- 
sion in these systems change the coherence and rms ampli- 
tude of the kHz QPOs? There is general agreement that the 
oscillation mechanism that produces the quasi-periodic vari- 
ability must be in the disk (see the references to theoretical 
models in iQ; it is very easy to come up with characteristic 
dynamical frequencies in the disk that match the observed 
frequencies of the kHz QPOs, whereas it appears more diffi- 
cult to have a "clock" somewhere else in the accretion flow. 
Nevertheless, from the early observations it became clear 
that emission from the disk alone cannot explain the rms 
amplitudes of the kHz QPOs, since in some cases the modu- 
lated luminosity in the QPOs is ~ 15% of the total emission, 
whereas at the same time the emission of the disk is ~ 10% 
or less. Furthermore, the steep increase of the rms ampli- 
tude with energy imply a large modulation o f the emitted 
flux, up to ~ 20% at ~ 25 — 30 keV (see e.g.. iBerger et alJ 
1996), at energies where the contribution of the disk is negli- 



gible. Hence, while the oscillation mechanism probably takes 
place in the disk, most likely the modulation mechanism is 
associated to the high-energy spectral component in these 
sources. 

The high-energy part of the X-ray spectrum of 
low-mass X-ray binaries ca n be described in terms 
of thermal Comptonization llSunvaev fc Titarchukl Il980l : 
White. Stella & Parmarl^l988). In low-luminosity Atoll 
sources I^H^nlgerfcTOm^^^rKlislll989^ . the hard spectral 
changes are explained as due to changes of the properties 
of this component as the average mass- accretion rate, M, 
chang es. For instance, for 4U 1608-52, lOierlihski fc Dona 
(2002) find that at low inferred M, the Comptonizing com- 
ponent has a relatively high temperature, T e > 20 keV, and 
is optically thin, r e ~ 1. The high-energy part of the spec- 
trum resembles a rather flat power law (power-law index 
~ 1.5 — 2), with a high-energy cut-off above ~ 100 keV 
jZhang et al.lll996h . At high inferred M. lQierliriski fc Done! 
( 2002) find that the Comptonizing component becomes 
cooler, with T e < 5 keV and the optical depth increases, 
T e >> 1. The high-energy part of the spectrum resembles 
a power law that is steeper (power-law index ~ 2 — 2.5) 
with a high-energy cut-off that is at lower energies than in 
the low M case. In the case o f high-luminosity Z-sources 
jHasinger fc van der Klislll989h . the Comptonizing compo- 
nent is relatively cool and optically thick as in the case of the 
low-l uminosity Atoll sources at high M dChristian fc Swankl 
Il997ft . 

Using a s i mple time-dependent Comptonization model, 
iLee fc Miller! |T998) calculated the spectrum of variability 
that would be produced by oscillations in the (i) injection 
rate of seed photons, (ii) density , and (iii) temperatur e of the 
Comptonizing medium (see also lStollman et alll987j) . They 
find that to reproduce the rms spectrum of the lower kHz 
QPO in 4U 1608-52, the variability must be mostly driven 
by an oscillation of the density of the Comptonizing medium. 
They also find that in the case of variations of the density of 
the Comptonizing medium, the largest rms variability in the 
1 — 10 keV range occurs as the optical depth of the medium 
is smallest (see their Figure 4, panel b). 

If the high-en ergy emission is d ue to Comptoniza- 
tion, the results of ILee fc Milled £l998), together with the 
global correlation between luminos ity and spectral hardness 
(Ivan Paradiis fc van der Klisl Il99i) . provide a possible ex- 
planation, at least qualitative, of the dependence of the rms 
amplitude of the kHz QPOs on luminosity and hardness. At 
low luminosity, corresponding to low mass accretion rate and 
hard spectra, the Comptonizing plasma is optically thin, and 
hence the amplitude of the variability is high. When mass 
accretion rate increases, the luminosity increases and the 
source becomes softer, and since the optical depth of the 
Comptonizing plasma increases, the amplitude of the vari- 
ability decreases. This would explain both the drop in rms 
amplitude of the kHz QPOs at high QPO frequency in indi- 
vidual sources, and that of the maximum rms amplitude in 
the sample of sources at high luminosity (low hardness). In 
individual sources, the rms amplitude of the kHz QPOs also 
decreases at low frequencies, which is more difficult to ex- 
plain in this simplified scenario. Also, in this scenario there 
is no straightforward explanation of the drop of the coher- 
ence of the lower kHz QPO at high frequencies in individual 
sources, and at high luminosities in the sample of sources. 
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iGilfanov. Revnivtsev. fc Molkovl J2003) have analyzed 
the rms energy spectra as a function of frequency of two low- 
mass X-ray binaries, the Z source GX 340+0 and the Atoll 
source 4U 1608-52. They find that the variability, including 
that in the frequency range of the kHz QPOs, is primarily 
due to variations of the luminosity of the boundary layer, 
with the emission of the accre tion d isk bein g much less vari - 
able. They also find (see also IGilfanov fc Revnivtsevll2005l) 
that in GX 340+0, the contribution of the boundary layer 
to the observed emission decreases as mass accretion rate 
increases, as the source moves along the Z-shaped track in 
the color-color diagr am from the Horizontal Branch to the 
Normal Branch (see lHasinger fc van der Klisi ri989. for an 
explanat ion of the branches in the color-color diagram of Z 
sources) . IGilfanov et al.l <l2003l) speculate that this decrease 
could be due either to obscuration of the boundary layer 
by a thickened accretion disk, a quantitative change in the 
structure of the boundary layer, or a complete disappearance 
of the boundary layer, when M ~ MEdd- In this scenario, 
the decrease of the rms variability of the kHz QPOs at high 
luminosities (high mass accretion rate) would be related to 
the lower relative contribution of the boundary layer, which 
produces the bulk of the QPO variability, to the total emis- 
sion. The obscuration mechanism cannot explain, however, 
the decrease of the coherence of the QPO at high luminos- 
ity. One way to explain the loss of coherence would be if the 
modulation mechanism in the boundary layer would damp 
out the oscillations. For instance, the flow of mass onto the 
boundary layer could be modulated more or less periodi- 
cally by the disc, and emission would proceed as a series 
of shots when mass reaches the boundary layer. Changes in 
the configuration o f the boundary layer, as those proposed 
bv IGilfanov et all l|2003f ). could increase the damping and 
reduce the lifetime of the QPO. This is of course just specu- 
lation; it remains to be seen whether it is possible to explain 
the decrease of QPO coherence by significant changes in the 
properties of the boundary layer. 

IB arret fc Verdenne I (1 19941) suggested that there is a 
critical luminosity for low-mass X-ray binaries with a 
neutron-star primary, L ~ 10 36 — 10 37 erg s ; sources be- 
low that level show a hard power-law component, whereas 
sources above that that level do not. Figure|^]shows a sort of 
dichotomy between the maximum rms amplitude and max- 
imum coherence of the kHz QPOs in the Z sources, for 
which L > 0.8 — 1.0I/Edd and the Atoll sources, for which 
L < 0.1 — 0.21/Edd- The separation, however, occurs at a lu- 
minosity that is ~ 5 — 10 times larger than suggested by 
iBarret fc Verdenne I lll994l) . Nevertheless, at least for the 
lower kHz QPO, the maximum coherence and maximum rms 
amplitude change significantly and rather smoothly within 
the Atoll class; this may indicate that, if the QPO proper- 
ties depend upon the emission of the power-law component, 
the relative importance of this compon ent may change in a 
less ab rupt manner than suggested bv IBarret fc Verdenne I 
Jl994h . 

The r max — r^ax correlation (Fig. [(Jl suggests that the 
same mechanism sets the amplitude of both QPOs. But in 
each source r max and r max occur at different luminosities 
(see Tables and |5J , and hence at different times, there- 
fore the mechanism cannot be acting simultaneously on both 
kHz QPOs. (In fact, in each source r max occurs at a higher 
frequency of the upper kHz QPO, v u than r max ; see, e.g., 



Ivan Straaten et al.l 1200 j . 120031) . If one ignores the case of 
4U 0614+09, the source with the lowest luminosity in the 
sample, also Q max is correlated with r max as wel l as with 
r max- From this it is tempting to infer that the mechanism 
that sets the rms amplitude of the kHz QPOs sets the coher- 
ence of the lower kHz QPO as well, at least for L>0.04Z/Edd- 
At low luminosities either this mechanism is not effective in 
setting Qmax, or another mechanism comes into play that 
counteracts the effect. The lack of a correlation between 
Qmax an d Qmaxi r max> ° r r max indicates that the mecha- 
nism that sets the latter three quantities does not have a 
significant effect on the coherence of the upper kHz QPO. 

As I showed in f3 for 4U 0614+09 either Q ma x is too 
low for r max , or r max is too high for Q^ ax . But since the 
maximum coherence for Le in 4U 0614+09 lies away from 
the J"max — Qmax as well as from the r max — Qmax correlations, 
it appears more likely that in this case Qf nax is too low both 
for the r max and r max , and not the other way around. This 
hypothesis is reinforced by the fact that a larger Q max would 
bring 4U 0614+09 closer to the other Atoll sources in Figures 
|^|and|SJ whereas lower r max and r max values would put 4U 
0614+09 in the gap of the correlations in those two Figures. 
If this is the case, it is unclear what mechanism could be 
responsible for the drop of coherence without a simultaneous 
drop of rms amplitude. 

Comparison of the upper right panels of Figures |3] and 
[7|suggests that the behaviors of Q u in individual sources (in 
this case 4U 1608-52) and Q max in the sample of sources 
as a function of hardness are also similar, in that both the 
coherence of the upper kHz QPO in individual sources, and 
the maximum coherence of the upper kHz QPO in the sam- 
ple of sources are more or less independent of the hard- 
ness of the source. However, the lack of correlation be- 
tween Qu and QPO frequency in 4U 1608-52 appears to 
be exceptional, and in other sources, both Atol l (see, e.g., 
Ivan Straaten et al.1 12002^: [Altamirano et alj|200 6 | ) as well as 
Z-sources fsee. e.g.^ ^ Tcter^Gis^tHnT^)97l : iHoman et alJ 
l2002Ujonker et al]l2000l.l2002l) . Q„ increases with QPO fre- 
quency, and hence decreases with hardness. From the data in 
Figure 0] it is not possible to discard a similar slow decrease 
of Q max vs - hardness. 

In black-hole systems, high-frequency QPOs have been 
observed in the range of tens up to ~ 450 Hz with typical 
rms amplitude and cohere nce values r ~ 1 — 3% and 
Q ~ 2 — 10, respectively |M prgan 1 Remillard, fc Greinerl 
1 19971: iRemillard et alJ Il999allbl 120021: ICui et alJ l200d: 
^^^nj™^^^^^^^^^^^^rriaB^e^|200j|_|^^^^^^3 
2001allbUHoman et al.ll200ll. l200l [20051: iMiller et al.ll200lf 
Klein- Wolt. Homan fc van der Klisl Eooi iBelloni et all 
200r^ ! The luminosities of these systems when high- 
frequency QPOs are observed are L ~ 1 — 3 x 10 38 
erg s" 1 (e.g in XTE J 1550-564; IHoman et ail 1200 ll : 
ISobczak etal]|200ll) or, for the typical mass of the black 
holes, Mbh ~ 10M©, L ~ 0.1Z/Edd- Compared to the 
neutron-star systems in Figure [3] in black-hole systems the 
rms amplitude and coherence of the high-frequency QPOs 
are too low for L/L-^dd of the source, but they would be 
consistent with the values for the Z sources at the same 
luminosity, i.e., the luminosities not being normalized to 
the source own Eddington luminosity (of course, to plot 
the black-hole sources in Figure [3] one would have to divide 
their luminosities by 2.5 x 10 38 erg s _1 , the Eddington 
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luminosity for a 1.9M© neutron star, as was done for the 
other sources in that plot). It remains to be seen whether 
in black-hole systems the modulation mechanism of the 
high-frequency QPOs would be similar to that of the kHz 
QPOs in neutron-star systems. 



5 CONCLUSION 

I study the maximum amplitude, r max , and maximum co- 
herence, Q m ax, of the kHz QPOs as a function of luminosity 
and hardness for a dozen low-mass X-ray binaries. I find 
that: 

(i) The maximum coherence of the lower kHz QPO, Qmax, 
first increases up to L ~ 0.04LEdd and then decreases with 
luminosity. 

(ii) The maximum coherence of the upper kHz QPO, 
Qmax, is independent of luminosity. 

(iii) The maximum rms amplitudes of both the lower and 
the upper kHz QPOs, r^ ax and r^ ax , respectively, decrease 
monotonically with luminosity 

(iv) Both r^ax and r^ax increase with the source hard- 
ness, Qmax first increases with hardness and then drops for 
the hardest source in the sample, and Q^ax is independent 
of hardness. 

(v) The relation between Q ma x and r max with luminosity 
in the sample of sources is similar to the relation between 
Q and r with QPO frequency in individual sources. The 
similarity extends also to hardness in the sample of sources 
and in at least one individual source, 4U 1608-52. 

(vi) The above argues against the interpretation that the 
drop of QPO coherence and QPO rms amplitude at high 
QPO frequency in individual sources is due to effects related 
to the innermost stable orbit around the neutron star in 
these systems. 

(vii) The drop of coherence and rms amplitude of the kHz 
QPOs, both in individual sources and in the sample could 
be produce by changes in the properties of the region in 
the accretion flow where the flux that produces the QPO is 
modulated. 
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